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On calculation of the neutron charge radius 
V.A. Karmanov^* 

'^Lebedev Physical Institute, Leninsky Prospekt 53, 119991 Moscow, Russia 

We show that the anomalous quark magnetic moments and relativistic effects in the 
nucleon wave function result in the correct value of the neutron charge radius. 

1. INTRODUCTION 

The neutron charge radius squared r|;„ (average of over the neutron charge distri- 
bution in coordinate space) is expressed through derivative of the neutron charge form 
factor Gp: 
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where F", are the form factors determining the neutron electromagnetic vertex: 
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Here mjv,u^ are the nucleon spinors, M is the nucleon mass and a^'^ = |(7^7^ — 1^1^)- 
According to (|1|), the neutron charge radius is represented as the sum of two contributions: 
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Hn = F^iS^) is the neutron magnetic moment. The contribution r^.^^^^ is the Foldy term [|l|, 
which appears due the generation of the electric field by the (anomalous) neutron magnetic 
moment because of its zitterbewegung. With the neutron magnetic moment /i„ = —1.91 



-0.126 /m^. The experimental value [@] 
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-0.113 ±0.005 fm^ is 



one finds: r^^i^y 

very close to r'^^^idy ^^"^ corresponds to r'l = 0.013 ± 0.005 fm'^, that is ten times smaller 
than Tp^i^y. 

Most of theoretical calculations, on the contrary, give r'l and r\^i^y to be close to each 
other in absolute values and of opposite signs, so, because of cancellation in the sum, 
is underestimated. In non-relativistic approximation, with the SU (6) symmetric wave 
function, the neutron charge form factor is identically zero (see, e.g., 0): G%{Q'^) = 0. 
To incorporate the relativistic effects, one usually proceeds from this wave function and 
takes into account the relativistic spins rotations by means of the Melosh matrices. They 
break the SU{Q) symmetry and give a non-zero value of r\^. However, as shown in 



* e-mails: karmanov@sci.lebedev.ru, karmanov@isn.in2p3.fr 



2 



the paper [|], the leading l/rn? correction, calculated in this way (m is the quark mass), 
results in = —r'^poidy = l/(3?7^^), so, in this approximation the value of r|;„ = rf + rp^i^^ 
is still zero. Similar result, by a different method, was found in Beyond the 
approximation, and rpoidy don't cancel each other, that can be seen from their analytical 
expressions up to the 

2 1 e 7(3' , _ 1 7e , _ 2e (3' 



where e = M — 3m, /5 is the parameter of the oscillator wave function (see eq. (0) below). 
We found by the methods of explicitly covariant light-front dynamics I^J^l- Since 
e ~ m and /S/m > 1, one can expect large high order corrections and invalidity of the 1/m 
decomposition at all. The numerical calculations without any approximation 
don't reveal a large effect and still underestimate r|;„: r\^{theor) ^ (0.65-=-0.70)r|;„(exp). 
In the paper [|l^ the value = —0.12 /m^ close to the experimental one was obtained, in 



the point form of dynamics, though there is a deviation of the nucleon magnetic moments. 

A reliable calculation of the neutron charge radius (as well as of the nucleon e.m. form 
factors) should incorporate many contributions, such as relativistic effects in the nucleon 
wave function, its nontrivial spin structure, exchange currents within nucleon, pionic 
fluctuations of constituent quarks [^, anomalous quark magnetic moments, etc. In this 
work we will show that the main discrepancy is easily removed, when the anomalous quark 
magnetic moments together with the relativistic effects are taken into account. The values 
of the quark magnetic moments are turned out to be close to those found in literature 
from other fits. Because of approximate and heuristic character of our consideration, this 
does not exclude other contributions, but shows importance of those which we take into 
account. 

2. SCALAR COUPLING AND ANOMALOUS QUARK MAGNETIC 
MOMENTS 

The full basis for the nucleon wave function contains sixteen spin structures 0. We will 
take the simplest one: 7(12,3) = (ui75f/cW2)(^3^Ar) = —7(21,3) corresponding to scalar 
coupling. Here Mi,2,3 are the (conjugated) quark spinors. The same coupling (but not 
only it), to construct the nucleon wave function, was used in the paper 0. 

We represent the symmetric spin-isospin part of the nucleon wave function in the form: 

ij = 1(12, 3)-e(12, 3) + /(23, l)-e(23, 1) + /(31, 2)i{31, 2), (5) 

where ^(12,3) = (^i^T2^2)(^3^Af) = —^21,3 is the isospin part with the zero pair isospin, 
^ is an isospinor. This wave function has the same non-relativist ic limit as the Melosh 
wave function, but it differs from the Melosh one in the relativistic domain. We assume 
that in neutron 123 = ddu. Calculating (^ for the particular isospinors C,d,d,u, we get: 
tpn = —I {23, 1) + /(31, 2). To simplify the further calculations, it is convenient to return 
to the canonical order 123. Using the Fierz identities, we transform tpn to the form: 

= -^(t'') 12 (7^75) 37V + ^{(^''ll2{cr^^-f5)3N, (6) 



^Note that the parities of u and UcU are opposite, hence ui'y^UcU2 is a true scalar, just due to 75. 
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where: {'y^)i2 = ui'y^UcU2, (7^75)3iv = u^lfilbUN, etc. Note that after Fierz transfor- 
mation the initial scalar coupling (|^) is transformed into the sum of vector and tensor 
couplings in (^. The total wave function is the product of the spin and the momentum 
parts: = ipn'^o- For ipQ we chose the oscillator form: 

i;o = Nexp[-My2P^), (7) 

where = {ki + k2 + k^Y is the invariant mass squared of 3g-system, A;i,2,3 are the 
on- mass shell quark four- momenta, N is found from the normalization condition 0. 

Let us introduce the anomalous quark magnetic moments k^, representing the quark 
e.m. vertex similarly to the nucleon one, eq. (^: 

U'qVUq = CqU'qYUq + l^'"' ^''^l' (^) 

We neglect the dependence of the quark form factors on Q^. The nucleon form factors are 
calculated in the explicitly covariant light-front dynamics , in the impulse approxima- 
tion, via plus-component of current, by a program of analytical calculation. The two-loop 
(five-dimensional) integrals are calculated numerically. 
In the 1/m? order the results can be obtained analytically: 

^ Qw? 2tm? 3m^' ^o'^f?/ 3,^2 gj^2 2,m? 
For the neutron charge radius we thus find: 

2 _ 2 I 2 _ 1 I ^ ?)K(i 

rEn-r^+rp^,,y- + 2,^2 + ^2 • l^Uj 

Note that when = = 0, we get r^^ = — So, with the wave function 
there is no cancellation in the 1 / order found in |^ with the Melosh wave function (we 
remind that both wave functions coincide with each other in non-relativistic limit). This 
cancellation is a peculiarity of the Melosh wave function. 

The fit of the nucleon static properties with m = 340 MeV, (3 = 700 MeV/c (that 
corresponds to (3"^ /m? ^ 4, < k >= 420 MeV/c), by numerical calculation not using the 
1/m decomposition, results in following anomalous quark magnetic moments: 

K„ = 0.13, Kd = -0.11. (11) 

With these values we get: 

fXn = -1.92, /ip = 3.09 and r|„ = -0.125 /m^ {rf = 0.002 /m^, r^^^^y = -0.127 /m^) 

The values of /i„ and r|;„ are rather close to the experimental ones, fip is by 10 % larger. 
With the wave function (|^), incorporating only one spin structure of sixteen, and with 
a simplest, oscillator parametrization we can hardly pretend to a higher accuracy. Note 
that the calculation by the approximate formula (0) gives rather close result: r|;„ = 
+ + ^) = -0.056 + (-0.045) ^ -0.1 frn^. From here one can see that half 
of r|;^ results from the relativistic effects in the nucleon wave function (the term — g^)- 
Another half results from the quark magnetic moments. 
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For comparison, we give Kg's found in other papers: Ku = 0.09, Kd = —0.15 @], k„ = 
-0.05 ^ 0.04, = 0.01 + -0.05 [|12[, k„ = 0.36, = -0.18 [gg. Though they differ 
from each other, these values (except for [|12|), as well as the values found in the present 
work, have a common features: {i) the absolute values of of the order of 0.1; {ii) 

is positive, Hd is negative. These features are enough to reproduce the experimental 
data, whereas the sensitivity of r^^^ to Kq is high. Putting Kg = or changing the sign 
of one or of both k^'s, we obtain the neutron magnetic moment and the charge radius 
incompatible with experiment (see Table |l|). 



Table 1 



Sensitivity of /i„ and to k^'s. The best fit: = +0.13, Kd = —0.11 





= +0.13 


Ku = Hu = —0.13 


K„ = +0.13 


K„ = —0.13 


Exp.: 


Kd = -0.11 


l^d = Kd = -0.11 


Kd = +0.11 


Kd = +0.11 


/in -1.91 


-1.92 


-1.46 -1.70 


-1.21 


-1.00 


rL,/m2 -0.113 


-0.125 


-0.07 -0.20 


+0.07 


-0.01 



We conclude that the anomalous quark magnetic moments (with the values close to 
already known in literature) and the relativistic effects in the nucleon wave function 
easily result in the correct neutron charge radius. 

I am grateful to L.Ya. Glozman for stimulating discussions attracting my attention to 
this problem. 
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^ In ref. ||T^, in contrast to eq. (^), the quark charge is separated as a factor at the front of both form 
factors. The original values in [y_3| are: k„ = = 0.54. Multiplying 0.54 by the quark charges 2/3 for u 
and —1/3 for d, we obtain the above values of Kg's from in the definition (||). 



